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Abstract 



^ ■ We study the properties of K and K mesons in nuclear matter at finite temper- 

C^^ ! ature from a chiral unitary approach in coupled channels which incorporates the s- 

and p- waves of the kaon-nucleon interaction. The in- medium solution accounts for 
^s) ■ Pauli blocking effects, mean-field binding on all the baryons involved, and tt and kaon 

r-^ , self-energies. We calculate K and K (off-shell) spectral functions and single particle 

C3 I properties. The K effective mass gets lowered by about —50 MeV in cold nuclear 

matter at saturation density and by half this reduction at T = 100 MeV. The p-wave 
contribution to the K optical potential, due to A, S and S* excitations, becomes sig- 
nificant for momenta larger than 200 MeV/c and reduces the attraction felt by the 
/\ • K in the nuclear medium. The K spectral function spreads over a wide range of en- 

ergies, reflecting the melting of the A(1405) resonance and the contribution of YN'"^ 
components at finite temperature. In the KN sector, we find that the low-density 
theorem is a good approximation for the K self-energy close to saturation density 
due to the absence of resonance-hole excitations. The K potential shows a moderate 
repulsive behavior, whereas the quasi-particle peak is considerably broadened with 
increasing density and temperature. We discuss the implications for the decay of the 
(j) meson at SIS/GSI energies as well as in the future FAIR/GSI project. 



PACS: 13.75.-n; 13.75. Jz; 14.20.Jn; 14.40.Aq; 21.65.+f; 25.80.Nv_ 

Keywords: Effective s-wave meson-baryon interaction, Coupled KN channels. Finite tem- 
perature. Spectral function, A(1405) in nuclear matter. 



1 Introduction 

The properties of hadrons and, in particular, of mesons with strangeness in dense matter 
have been a matter of intense investigation over the last years, in connection to the study 
of exotic atoms [1] and the analysis of heavy-ion collisions (HIC's) [2]. 

At zero temperature, the study of the K interaction in nuclei has revealed some interest- 
ing characteristics. First, the presence below the KN threshold of the A(1405) resonance 
gives rise to the failure of the Tp approximation for the K self-energy. Whereas the KN 
interaction is repulsive at threshold, the phenomenology of kaonic atoms requires an at- 
tractive potential. The consideration of Pauli blocking on the intermediate KN states [3,4] 
was found to shift the excitation energy of the A (1405) to higher energies, hence chang- 
ing the real part of the KN amplitude from repulsive in free space to attractive in a 
nuclear medium already at very low densities. Further steps were taken to account for self- 
consistency in the evaluation of the K self-energy [5] as well as for relevant medium effects 
on the intermediate meson-baryon states [6] , which results in a moderate final size of the 
attractive potential as well as in a sizable imaginary part associated to several in-medium 
decay mechanisms [6-9]. Several studies of the K potential based on phenomenology of 
kaonic atoms have pointed towards a different class of deeply attractive potentials [10]. 
Unfortunately, the present experimental knowledge is unable to solve this controversy, as 
both kinds of potentials fairly describe the data from kaonic atoms [10, 11]. 

A different direction in the study of this problem was given in [12], where a highly 
attractive antikaon-nucleus potential was constructed leading to the prediction of narrow 
strongly bound states in few body systems [12-14]. This potential has been critically dis- 
cussed in [15] because of the omission of the direct coupling of the vrS channel to itself, the 
assumption of the nominal A(1405) as a single bound K state, the lack of self-consistency in 
the calculations and the seemingly too large nuclear densities obtained of around ten times 
normal nuclear matter density at the center of the nucleus. Experiments devised to the 
observation of deeply bound kaonic states, measuring particles emitted after the absorption 
of K~ in several nuclei, reported signals that could actually be interpreted from conven- 
tional nuclear physics processes. The experimental observations could be explained simply 
either by the two-body absorption mechanism raised in [15], without [16-19] or with [20,21] 
final state interactions, or coming from a three-body absorption process [22-24] . Actually, 
recent improved few-body calculations using realistic KN interactions and short-range cor- 
relations [25-28] predict few-nucleon kaonic states bound only by 50-80 MeV and having 
large widths of the order on 100 MeV, thereby disclaiming the findings of Refs. [12-14]. 

Relativistic heavy-ion experiments at beam energies below 2AGeV [29,30] is another 
experimental scenario that has been testing the properties of strange mesons not only in a 
dense but also in a hot medium. Some interesting conclusions have been drawn compar- 
ing the different theoretical transport-model predictions and the experimental outcome on 
production cross sections, and energy and polar-angle distributions [29]. For example, de- 
spite the significantly different thresholds in binary A^A^ collisions, there is a clear coupling 
between the K~ and K^ yields since the K~ is predominantly produced via strangeness 
exchange from hyperons which, on the other hand, are created together with K^ mesons. 



Furthermore, the K^ and K^ mesons exhibit different freeze-out conditions as the K~ 
are continuously produced and reabsorbed, leaving the reaction zone much later than the 
K^ mesons. However, there is still not a consensus on the influence of the kaon-nucleus 
potential on those observables [31]. 

Compared to the KN interaction, the KN system has received comparatively less at- 
tention. Because of the lack of resonant states in the S* = +1 sector, the single-particle 
potential of kaons has usually been calculated in a Tp approximation, with a repulsion of 
around 30 MeV for nuclear matter density (the information on T is taken from scattering 
lengths and energy dependence is ignored). However, a recent analysis of the KN interac- 
tion in the Jiilich model has demonstrated that the self-consistency induces a significant 
difference in the optical potential with respect to the low-density approximation at satu- 
ration density, and the kaon potential exhibits a non-trivial momentum dependence [32]. 

A precise knowledge of kaon properties in a hot and dense medium is also an essential 
ingredient to study the fate of the meson. Electromagnetic decays of vector mesons offer 
a unique probe of high density regions in nuclear production experiments and HIC's [2]. 
The meson predominantly decays into KK, which are produced practically at rest in the 
center of mass frame, each carrying approximately half of the mass of the vector meson. 
Such a system is highly sensitive to the available phase space, so that small changes in 
the kaon effective masses or the opening of alternative baryon-related decay channels may 
have a strong repercussion in the 0-meson mass and decay width. The analysis of the mass 
spectrum of the decay products in dedicated experiments has drawn inconclusive results 
since the long-lived vector meson mostly decays out of the hot/dense system [33-38]. In 
addition, despite the sizable modifications predicted in most theoretical studies [39-46], the 
current experimental resolution typically dominates the observed spectrum from decays. 
Still, recent measurements of the transparency ratio in the nuclear photoproduction 
reaction by the LEPS Collaboration [36] have shed some light on the problem and seem 
to indicate an important renormalization of the absorptive part of the nuclear potential, 
as it was suggested in the theoretical analysis of Ref. [47]. The unprecedented precision 
achieved by the CERN NA60 Collaboration in the analysis of dimuon spectrum data from 
In-In collisions at 158 AGeV [48] , as well as the advent of future studies of vector meson 
spectral functions to be carried out at the HADES [49] and CBM [50] experiments in the 
future FAIR facility, advise an extension of our current theoretical knowledge of the 
spectral function to the {hb,T) plane, and hence, of the K and K properties at finite 
temperature and baryon density. 

In this work we evaluate the K and K self-energy, spectral function and nuclear optical 
potentials in a nuclear medium at finite temperature. We follow the lines of Refs. [6,9,51,52] 
and build the s-wave KN and KN T-matrix in a coupled channel chiral unitary approach. 
Medium effects are incorporated by modifying the intermediate meson-baryon states. We 
account for Pauli blocking on intermediate nucleons, baryonic binding potentials and me- 
son self-energies for pions and kaons. The latter demands a self-consistent solution of the 
K and K self-energies as one sums the kaon nucleon scattering amplitude over the occu- 
pied states of the system, whereas the T-matrix itself incorporates the information of the 
kaon self-energies in the intermediate meson-baryon Green's functions. The interaction 



in p-wave is also accounted for in the form of YN~'^ excitations, which lead to a sizable 
energy dependence of the K self-energy below the quasi-particle peak. Finite temperature 
calculations have been done in the Imaginary Time Formalism in order to keep the re- 
quired analytical properties of retarded Green's functions which, together with the use of 
relativistic dispersion relations for baryons (and, of course, for mesons), improves on some 
approximations typically used in former works. The organization of the present article goes 
as follows: in Sect. [2] we develop the formalism and ingredients on which the calculation 
is based. Sect. |3]is devoted to the presentation of the results. The K and K self-energies 
and spectral functions are discussed in Sects. [3TT] and 13. 2[ respectively. We devote Sect. 13.31 
to the discussion of momentum, density and temperature dependence of K and K nuclear 
optical potentials. Finally, in Sect. |4]we draw our conclusions as well as the implications 
of the in-medium properties of kaons at finite temperature in transport calculations and 
meson phenomenology. We also give final remarks pertaining to the present and future 
works. 

2 Kaon nucleon scattering in hot nuclear matter 

In this section we discuss the evaluation of the effective kaon nucleon scattering amplitude 
in a dense nuclear medium, extending the unitarized chiral model for K of Refs. [6, 9] 
to account for finite temperature. This allows us to obtain the in-medium K and K 
self-energy, spectral function and nuclear optical potential. We follow closely the lines of 
Ref. [51], where a similar study was reported for open-charm mesons in hot nuclear matter. 

2.1 s-wave kaon nucleon scattering and kaon self-energy 

The kaon nucleon interaction at low energies has been successfully described in Chiral 
Perturbation Theory (xPT) [53-57]. The lowest order chiral Lagrangian which couples the 
octet of light pseudoscalar mesons to the octet of 1/2+ baryons is given by 

/:i^) = (BiYV^B) - M{BB) 

+^D {bYi, {u,, B}) + ^F (57^^75 K, B]) , (1) 

where B is the 5'[/(3) matrix for baryons, M is the baryon mass, u contains the $ matrix 
of mesons and the symbol ( ) denotes the flavour trace. The S'f/(3) matrices appearing in 
Eq. ([1]) are standard in notation and can be found, for instance, in [52]. The axial- vector 
coupling constants have been determined in [58] and read D = 0.8 and F = 0.46. 

Keeping at the level of two meson fields, the covariant derivative term in Eq. ([T]) provides 
the following interaction Lagrangian in s-wave, 

4^) = /5z7^^[(<l>9^<l>-9^$<I>)fi-fi($a,,$-9^$$)]\ . (2) 



Expanding the baryon spinors and vertices in M^^ the following expression of the s-wave 
kaon nucleon tree level amplitude can be derived 

with Mb^ and Ei the mass and energy of the baryon in the ith channel, respectively. The 
coefficients Cij form a symmetric matrix and are given in [52]. The meson decay constant 
/ in the s-wave amplitude is taken as / = I.IS/tt. This renormalized value provides a 
satisfactory description of experimental low energy KN scattering observables (such as 
cross sections and the properties of the A(1405) resonance) using only the interaction from 
the lowest order chiral Lagrangian plus unitarity in coupled channels. We have considered 
the following channels in the calculation of the scattering amplitude: in the strangeness 
sector S* = — 1 we have KN, vrS, rjA and KS for isospin / = 0; KN, ttA, ttS, t^S and KS 
for 1 = 1. For S* = 1, there is only one single channel, KN, for each isospin. 

In Refs. [6, 9, 52] unitarization of the tree level amplitudes in coupled channels was 
implemented, which extends the applicability of xPT to higher energies and in particular 
allows to account for dynamically generated resonances. In particular, the A(1405) shows 
up in the unitarized s-wave KN amplitude. Following [52], the effective kaon-nucleon 
scattering amplitude is obtained by solving the Bethe-Salpeter equation in coupled channels 
(in matrix notation), 

T = V + VGT, (4) 

where we use the s-wave tree level amplitudes as the potential (kernel) of the equation, 
Vfj, and 

r(n = - f -^ -^ ^ 1 ..^ 

^^^'^ ' J (2 ny E,{-q) ^s-q,- E,{-q) + ie ql - q^ - mj + is ^ ' 

stands for the intermediate two-particle meson-baryon Green's function of channel i {G 
is diagonal). In principle, both V and T enter off-shell in the momentum integration 
{VGT term) of the meson-baryon loop. However, as it was shown in Refs. [9,52], the 
(divergent) off-shell contributions of V and T in the s-wave interaction can be reabsorbed 
in a renormalization of the bare coupling constants and masses order by order. Therefore, 
both V and T can be factorized on-shell out of the meson-baryon loop, leaving the four- 
momentum integration only in the two-particle meson-baryon propagators. An alternative 
justification of solving the Bethe-Salpether equation with on-shell amplitudes may be found 
in the framework of the N/D method, applied for meson-meson interactions in Ref. [59] and 
for meson-baryon interactions in Ref. [60]. We are thus left with a set of linear algebraic 
equations with trivial solution, 

T = [1 - VG]-W . (6) 

The meson baryon loop function, Gi, needs to be regularized. We apply a cut-off in the 
three momentum of the intermediate particles, which provides a simple and transparent 
regularization method for in- medium calculations, cf. [6,9]. 



In order to obtain the effective s-wave K{K)N amplitude in hot and dense matter, we 
incorporate in the loop functions the modifications on the properties of the mesons and 
baryons induced by temperature and density. 

In the Imaginary Time Formalism (ITF), the baryon propagator in a hot medium is 
given by: 

where ium = i(2m + l)7rT + fis is the fermionic Matsubara frequency, with /i^ the baryon 
chemical potential, and Eb is the baryon single particle energy, which, in the case of nucle- 
ons and singly strangeness hyperons, will also contain the medium binding effects obtained 
within a temperature dependent Walecka-type a — u model (see Ref. [61]). According to 
this model, the nucleon energy spectrum in mean-field approximation is obtained from 



Ej,ip, T) = ^p + M*^{Ty + S'' , (8) 

with the vector potential S^ and the effective mass M'^{T) given by 

M*^{T) = Mn- S^ with S^ = (^\ ps , (9) 

\ms/ 

where m^ and m^ are the meson masses {rris = 440 MeV, m^ = 782 MeV), while Qs and 
Qv are the scalar and vector density dependent coupling constants. These constants are 
obtained by reproducing the energy per particle of symmetric nuclear matter at T = 
coming from a Dirac-Brueckner-Hartree-Fock calculation (see Table 10.9 of Ref. [62]). The 
vector (p) and scalar (ps) densities are obtained by momentum integration, namely 

m=4/^nS^Hp,T), (10) 

of the corresponding vector [nAr(p, T)] and scalar [n%{'p,T)] density distributions, which 
are defined in terms of the nucleon Fermi-Dirac function as 

nN{p,T) = — — -- (11) 

1 + exp [{En{p, T) - Pb)/T\ 

and 

respectively. The quantities En{'p,T),M'^(T) and pb are obtained simultaneously and 
self-consistently for given p and T and for the corresponding values of Qs and g^. 
The hyperon masses and energy spectra. 



Ey{p,T) = dp + M^{TY + S^ , (13) 



can be easily inferred from those for the nucleon as 

6 xnivJ o 

M^{T) = My-^y = My-- (—) ps 



3 Vm, 

= My-'^iMr.-MUT)) . (14) 

Here we have assumed that the a and uj fields only couple to the u and d quarks, as in 
Refs. [63,64], so the scalar and vector coupling constants for hyperons and charmed baryons 
are 

9v = -^dv, gj = g^s- (15) 

In this way, the potential for hyperons follows the simple light quark counting rule as 
compared with the nucleon potential: Vy = 2/3 Vn- As reference, we quote in Table 
[1] the nucleon and hyperon single particle properties for three densities (0.25po5 Po and 
2po), where po = 0.17 fm~^ is the normal nuclear matter saturation density, and four 
temperatures (T = 0, 50, 100 and 150 MeV). The hyperon attraction at p = po and T = 
MeV is about —50 MeV, the size of which gets reduced as temperature increases turning 
even into repulsion, especially at higher densities. This behavior results from the fact 
that the temperature independent vector potential takes over the strongly temperature- 
dependent scalar potential which decreases with temperature. We note that the quark- 
meson coupling (QMC) calculations of Refs. [63,64], performed at T = 0, obtained a 
somewhat smaller scalar potential (about half the present one) for the A and E baryons 
due to the inclusion of non-linear terms associated to quark dynamics. To the best of our 
knowledge, no temperature effects have been studied within this framework. 
The meson propagator in a hot medium is given by 

Dm(o;„, g; T) = , (16) 

where iti;„ = Ylm^T is the bosonic Matsubara frequence and IiMi!jJni<i]T) is the meson 
self-energy. Note that throughout this work we set the mesonic chemical potential to zero, 
since we are dealing with an isospin symmetric nuclear medium with zero strangeness. We 
will consider in this work the dressing of pions and kaons. 

An evaluation of the in-medium self-energy for pions at finite temperature and baryonic 
density was given in the Appendix of Ref. [8], which generalized the zero temperature 
evaluation of the pion self-energy from Refs. [65, 66] by incorporating thermal effects. We 
recall that the pion self-energy in nuclear matter at T = is strongly dominated by the 
p-wave coupling to particle-hole {ph) and A-hole (A/i) components (a small, repulsive s- 
wave contribution takes over at small momenta), as well as to 2p-2h excitations, which 
account for pion absorption processes, and short range correlations. We come back to the 
pion self-energy in Sec. 12. 3[ where we improve on some approximations in previous works. 



Table 1: a — u model at finite temperature 

p[fm~^] T[MeV] /ii?[MeV] M*^[MeV] S^[MeV] MX[MeV] SX[MeV] M^[MeV] E|;[MeV] 

781 121 1010 81 1088 81 

793 121 1019 81 1096 81 

805 121 1026 81 1104 81 

815 121 1033 81 1110 81 

579 282 872 188 950 188 

605 282 893 188 970 188 

634 282 913 188 990 188 

659 282 929 188 1006 188 

443 422 787 281 865 281 

470 422 803 281 881 281 

510 422 830 281 907 281 

545 422 853 281 930 281 

In the case of the kaons, the self-energy receives contributions of comparable size from 
both s- and p-wave interactions with the baryons in the medium. We evaluate the s-wave 
self-energy from the effective K{K)N scattering amplitude in the medium, a procedure 
which, as will be shown explicitly in the following, must be carried out self-consistently. 
The p-wave part of the kaon self-energy will be discussed separately in the next section. 

The evaluation of the effective K{K)N s-wave scattering amplitude in the hot medium 
proceeds by first obtaining the meson-baryon two-particle propagator function (meson- 
baryon loop) at finite temperature in the ITF, Qmb- Given the analytical structure oiQMB 
it is convenient to use the spectral (Lehmann) representation for the meson propagator. 
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DM{iOn,q;T) = duj 



SM{,u;,q;T) 



iu„ — u 



iuJn — id Jo iUn + UJ 

where Sm, S^ stand for the spectral functions of the meson and its corresponding anti- 
particle. The separation in the second line of Eq. ( TT7I) reflects the retarded character of the 
meson self-energy and propagator, Imn(— go, (f; T) = — Imn(go, ^; 7"). Due to strangeness 
conservation, K and K experience markedly different interactions in a nuclear medium [61 , 
which justifies using a different notation {Sm, S^) for the spectral functions in Eq. (fTTll ^l. 
Combining Eqs. (iT6l) and (IT71) . conveniently continued analytically from the Matsubara 
frequencies onto the real energy axis, one can write 

SM{^,q]T) = Im DM{uj,q;T) = ^ TTl ^^FVJ^ ' ^^^^ 



^In the case of pions, for instance, in an isospin symmetric nuclear medium, all the members of the 
isospin triplet (tt^, 7r°) acquire the same self-energy and one can write 5'^(— (?0; T, T) — —S-jriqa, T, T) with 
the subsequent simplification of Eq. p?|) . 
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Applying the finite-temperature Feynman rules, the meson-baryon loop function in the 
ITF reads 



gMB{W^,P;T) 



-T 



cP'q 



E 



1 



X / duo 
h 



(27r)3 -^ IWm - lOOn - Eb{P - q, T) 
SM{uJ,q;T) SM{uj,q;T)\ 



lUJr. 



IjJ 



\uJn + ^0" 



(19) 



where P is the external total three-momentum and Wm an external fermionic frequency, 
iW^m = i(2m+l)7rT+/iB. Note that we follow a quasi-relativistic description of baryon fields 
all throughout this work: in Eq. flT9|) the negative energy part of the baryon propagator 
has been neglected but we use relativistic dispersion relations. All the Dirac structure 
is included in the definition of the tree level amplitudes. The Matsubara sums can be 
performed using standard complex analysis techniques for each of the two terms in the 
meson propagator and one finds 



Qmb{W^,P-T) 



Sq 

(2^ 



dio 



SM{uj,q]T) 



l-nB{P-q,T) + f{u;,T) 
iW^-uj-EB{P-q,T) 



+ SM{uj,q;T)- 



nBiP-q,T) + f{u;,T) 
iWm + co-EBiP-q,T) 



(20) 



with f{u),T) = [exp(co'/T) — 1] ^ the meson Bose distribution function at temperature T. 
The former expression can be analytically continued onto the real energy axis, Gmb{Po + 
ie ,P;T) = QMB{iWm — > Po + ^£ ,P',T), cf. Eq. ( !20|) . With these medium modifications the 
meson-baryon retarded propagator at finite temperature (and density) reads 
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(21) 



for K{K)N states and 

G^Y{Po + ie,P;T) 



d^q 



M 



Y 



+ 



i^Tif Ey{P - q,T) 

X) 

duST,{uj, q,T) 



l + fico,T) 



Po + ie-u-EYiP-q,T) 



Po + ie + u-EYiP-q,T) 
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f22) 



for ttA or ttS states, where P = {Po,P) is the total two-particle niomentum and g* is 
the meson three-momentum in the nuclear medium rest frame. Note that for consistency 
with the free space meson-baryon propagator, given in Eq. ([5]), we have included the 
normalization factor Mb/ Eb in the baryon propagator. We have explicitly written the 
temperature dependence of the baryon energies indicating that we account for mean-field 
binding potentials as discussed above. The second term in the KN loop function typically 
provides a small, real contribution for the energy range in Pq we are interested in. In 
order to simplify the numerical evaluation of the self-consistent KN amplitude we replace 
Sk{^,(1',T) by a free-space delta function in Eq. ( 12T1) . This approximation is sensible as 
long as the K spectral function in the medium still peaks at the quasi-particle energy 
and the latter does not differ much from the energy in vacuum, as we will confirm in 
Sect. 131 In Eq. (12T]) we have neglected the kaon distribution function, since we expect Bose 
enhancement to be relevant only for the lightest meson species in the range of temperatures 
explored in the present study, T = - 150 MeV. 

The 7cY loop function, in particular, incorporates the 1 + f{uJ,T) enhancement factor 
which accounts for the contribution from thermal pions at finite temperature, cf. Eq. fl22l) . 
In this case, we have neglected the fermion distribution for the participating hyperons, 
which is a reasonable approximation for the range of temperature and baryonic chemical 
potential that we have studied (cf. Table [1]). 

In the case of r]A, vjTj and K'E intermediate states, we simply consider the meson 
propagator in free space and include only the effective baryon energies modified by the 
mean-field binding potential, namely 

GKPo + i^, P-,T)= lf\. \^, J^' ^ . . (23) 

J (27:)^ 2uJi{q)Ei{P-q,T) Po + ie-iJi{q)-Ei{P-q,T) 

The latter channels are less relevant in the unitarization procedure of the s-wave scattering 
amplitude. They are important to maintain SU(3) symmetry through using a complete ba- 
sis of states in the coupled-channels procedure, as well as for producing a better description 
of branching ratios between the various scattering transitions at threshold. However, the 
width and position of the A (1405) are basically determined from the unitarized coupling of 
^A^ and vrS channels [52]. In addition, the changes that kaons [32,67] and t] mesons [68,69] 
experience in the medium at moderate densities are comparably weaker than for vr and K, 
which justifies the simplification adopted here. 

As mentioned above, all meson-baryon loop functions in our approach are regularized 
with a cut-off in the three-momentum integration, gmax- We adopt here the regularization 
scale that was set in [52] , where the KN scattering amplitude was evaluated in free space, 
leading to a remarkable description of several KN scattering observables and the dynamical 
generation of the A(1405) resonance with a single parameter, g^ax = 630 MeV/c. 

Finally, we obtain the in-medium s-wave K{K) self-energy by integrating Tj^,^j^-^^ over 
the nucleon Fermi distribution at a given temperature, 

n^(x)(?o, q; ^) = / (0j ^Nip^ T) [4'^^^(Po, P; t) + ^T'i=^]^{P,, P- T)] , (24) 

10 



where Pq = qo + En{p, T) and P = q + p are the total energy and momentum of the 
K{K)N pair in the nuclear medium rest frame, and q stands for the momentum of the 
K{K) meson also in this frame. The kaon self-energy must be determined self-consistently 
since it is obtained from the in-medium amplitude, Tj^(k)n^ which requires the evaluation 
of the K[K)N loop function, Gif(x)Ar) and the latter itself is a function of n/;-(^)(go, q] T) 
through the kaon spectral function, cf. Eqs. flTSj) . fl2Tl) . Note that Eq. fl24l) is valid in cold 
nuclear matter. In the Appendix we provide a derivation of the s-wave self-energy from 
the kaon nucleon T-matrix in ITF and give arguments for the validity of Eq. fl2^ as a 
sensible approximation. 

2.2 j9-wave kaon self-energy 

The main contribution to the p-wave kaon self-energy comes from the A and S pole terms, 
which are obtained from the axial-vector couplings in the Lagrangian {D and F terms in 
Eq. ([T])). The S*(1385) pole term is also included explicitly with couplings to the kaon- 
nucleon states which were evaluated from SU(6) symmetry arguments in [43]. 

In Ref. [9] , the combined s- and p-wave K self-energy was obtained self-consistently in 
cold nuclear matter by unitarizing the tree level amplitudes in both channels. Unitarization 
of the p-wave channel, however, did not provide dramatic effects over the tree level pole 
terms. Those were used in [43] to evaluate the p-wave self-energy. The latter is built 
up from hyperon-hole (Yh) excitations and it provides sizable strength in the K spectral 
function below the quasi-particle peak and a moderate repulsion (as it is expected from 
the excitation of subthreshold resonances) at nuclear matter density. Medium effects and 
unitarization actually move the position of the baryon pole with respect to that in free 
space [9]. We incorporate this behaviour here in an effective way through the use of 
baryon mean field potentials and re-evaluate the relevant many-body diagrams for the 
p-wave self-energy at finite temperature and density, which considerably simplifies the 
numerical task. We shall work in the ITF, improving on some approximations of former 
evaluations [8,43,44]. Hence, the total K and K self-energy will consist of the sum of the 
s- and p-wave contributions described here and in the former section. Note that the p-wave 
self-energy enters the self-consistent calculation of the s-wave self-energy as the K and K 
meson propagators in the intermediate states on the s-wave kaon-nucleon amplitude are 
dressed with the total self-energy at each iteration. 

We can write the K p-wave self-energy as the sum of each of the Yh contributions, 

n^(?o, q;T) = - V^^^ q^ fl{qo, q) Uan-^ (go, q; T) 
3 ~ 

+ 2^RNj:*^^M*i'lo,?)U^*N~^{qo,q;T) , (25) 



where Uy stands for the YN ^ Lindhard function at finite temperature and baryonic 

KNY 



density, and V^^ry represents the KNY coupling from the chiral Lagrangian in a non- 
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relativistic approximation (leading order in a M^^ expansion) and includes the required 
isospin multiplicity. They can be found, for instance, in Ref. [43]. The fy factors account 
for relativistic recoil corrections to the KNY vertices, which improve on the lowest order 
approximation and still allow to write the self-energy in a simple form, where all the 
dynamical information from the p-wave coupling is factorized out of the momentum sum 
in the fermionic loop. These factors read 

/ls(?o,<?) = [2MA,s + 2go(MA,s-M^)-g2 + 2M,v^A,s(g)]/4Mjv^A,s(g') , 
/|.(go,g) = (1-go/MsO' ■ (26) 

We have also accounted for the finite size of the vertices by incorporating phenomenological 
hadronic form factors of dipole type via the replacement q'^ -^ Fxiq"^) q"^ with F^(g^) = 
(A|^/(A|; + q^)y, where A^ = 1050 MeV. We provide exphcit expressions for Uyn-'^ in 
Appendix [Al 

In Refs. [43,44] the K and K self-energy in p-wave was obtained for cold nuclear mat- 
ter, and it was extended to the finite temperature case in [8]. We would like to present 
here a comparison of our results for different approximations as the ones used in former 
evaluations, particularly in the T —>■ limit. In Fig. [1] we present the imaginary part of the 
K p-waye self-energy as a function of the kaon energy, evaluated at nuclear matter density 
and two different kaon momenta. In the upper panels the dashed lines have been obtained 
by using, in Eq. fl2^ . the standard non-relativistic evaluation of the YN~^ Lindhard func- 
tion [43,44], the analytical expression of which is also quoted in Appendix |Al The solid 
lines correspond to our calculation in Eq. (133!) for T = MeV. Both results agree quite well 
which ensures that our calculation has the correct T — ;> limit. The observed differences, 
particularly the threshold energies at which each YN~^ component is open/closed, are 
related to the use of non-relativistic baryon energies (in dashed lines) and some further 
approximations that we discuss below. 

At T = 100 MeV (lower panels), the dashed lines correspond to the evaluation of the 
Lindhard function in [8], which extends the zero-temperature, non-relativistic calculation 
by replacing the nucleon occupation number, 6{pf — p), by the corresponding Fermi-Dirac 
distribution, n]^{p,T). Note that, in our result (solid lines), the Matsubara sum automat- 
ically generates a non-vanishing term proportional to the hyperon distribution function, 
cf. Eqs. (I31II33P in Appendix [Al Despite the absence of a crossed kinematics mechanism (as 
for instance in ph and A/i excitations for pions) this guarantees that the imaginary part 
of the K self-energy identically vanishes at zero energy (go = 0) as it follows from the (re- 
tarded) crossing property of the thermal self-energy, Imn^(— go, g; 7") = — Inin^(go, g; 7") 
Q The two calculations at T = 100 MeV exhibit some non-trivial differences. In the 
finite-T extension of the non-relativistic result the YN~^ structures and thresholds are 
more diluted and the S component is not resolved even at low momentum, whereas in the 
relativistic result the three components are clearly identifiable at g = 150 MeV/c. In both 
calculations the strength of the imaginary part extends to lower energies so that the energy 



^At zero energy (go = 0) the K and K modes cannot be distinguished and thus lnin^(0, g;r) = 
— Imn^(0,g;T) = 0. This is accomphshed exactly in Eq. p3|) in Appendix \K\ for go = 0. 
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Figure 1: Imaginary part of the K p-wave self-energy from different evaluations of the 
YN~^ Lindhard function (see text) at p = po- The upper panels correspond to T = 0, 
whereas in the lower panels T = 100 MeV. 



gap in cold nuclear matter is absent here. However, as mentioned before, the relativistic 
calculation in the ITF vanishes exactly at go = whereas the non-relativistic, finite-T 
extended result does not. In addition, the standard non-relativistic calculations (at T = 
and finite T) [8,43,44] omit p^ terms proportional to {My^ — M^^) in the energy balance 
of the two-particle propagator, which are responsible for cutting the available YN~^ phase 
space for high external energies (go)- As a consequence, the dashed lines in the lower pan- 
els exhibit a high energy tail which is attenuated by the nucleon distribution whereas the 
relativistic result has a clear (temperature dependent) end point. This happens for each 
of the excited hyperon components independently and it is particularly visible for the S*. 
The A and S tails mix with each other and are responsible for the washing out of the S 
structure at finite momentum. 

Finally, the K p-wave self-energy can be obtained following the same procedure as 
above. The excitation mechanisms in this case correspond exactly to the crossed kinematics 
of those of the K. Taking advantage of the crossing property of the thermal self-energy, 
we obtain the K self-energy from the K one replacing go — ^ —go in 11^ (go, g; T) (modulo 
a sign flip in the imaginary part). The crossed kinematics causes the YN~^ excitations 
to be far off-shell. In cold nuclear matter, for instance, 11^ is real and mildly attractive. 
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At finite temperature, though, the fermion distributions of the nucleon and hyperon can 
accommodate low-energy (off-shell) kaons and 11^ receives a small imaginary part which 
rapidly decays with increasing energy. 

2.3 Pion self-energy 

We briefly discuss here the relevant many-body mechanisms that modify the pion propa- 
gator, which enters the evaluation of the in- medium KN amplitude, cf. Eqs. ([6]), fl22|) . In 
cold nuclear matter, the pion spectral function exhibits a mixture of the pion quasi-particle 
mode and p/i, A/i excitations [65]. The meson-baryon chiral Lagrangian in Eq. ([T]) pro- 
vides the ttNN p-wave vertex, while the nNA vertex can be determined from the standard 
non-relativistic derivation of the Raritta-Schwinger interaction Lagrangian. However, we 
shall use phenomenological irNN and vriVA coupling constants determined from analysis 
of pion nucleon and pion nucleus reactions. Their values are /Ar/m,r = 0.007244 MeV~^ 
and /a/Ztv = 2.13. The lowest order p-wave pion self-energy due to ph and Ah excitations 
then reads 




n^iviv-i+^Aiv-i(?o, q;T) = [ — q pMN-^igo, q; t) + UAN~^{qo, g; t)] , (27) 



where the finite temperature Lindhard functions for the ph and Ah excitations are given 
in Appendix \M Note that, for convenience, we have absorbed the nNA coupling in the 
definition of U/^^n-i. 

The strength of the considered collective modes is modified by repulsive, spin-isospin 
A^A^ and A^A short range correlations [70], which we include in a phenomenological way 
with a single Landau- Migdal interaction parameter, g' = 0.7. The RPA-summed pion 
self-energy then reads 

^,, .^, {{^fFA?') q' [UNN-^jqo, g; T) + [/^^-^(go, g; T)] 

ii^igo,g;jj- 7t\2 > i-^^J 

1- (^) g' [f/iV7v-i(go,g1T) + f/A7v-i(go,g;T)] 

which also contains the effect of the same monopole form factor at each nNN and ttA^A 
vertex as used in T = studies, namely Fj.^q'^) = (A^ — m^)/[A^ — (go)^ + g^], with 
A^ = 1200 MeV, as is needed in the empirical study of A^A^ interactions. Finally, for 
consistence with former evaluations of the pion self-energy, we have also accounted for 
one-body s-wave scattering and 2p2h mechanisms, following the results in Refs. [66,71,72], 
which we have kept to be the same as in T = 0. 

In Fig. [2] we show the pion spectral function at normal nuclear matter density for two 
different momenta. At T = (upper panels) one can easily distinguish the different modes 
populating the spectral function. At low momentum, the pion quasi-particle peak carries 
most of the strength together with the ph structure at lower energies. The Ah mode starts 
to manifest to the right hand side of the pion quasi-particle peak. Note that the pion 
mode feels a sizable attraction with respect to that in free space. At higher momentum, 
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the excitation of the A is clearly visible and provides a considerable amount of strength 
which mixes with the pion mode. As a consequence, the latter broadens considerably. The 
solid lines include also the contributions from the pion s-wave self-energy and two-body 
absorption. These mechanisms, especially the latter, generate a background of strength 
which further broadens the spectral function, softening in particular the pion peak at low 
momentum and the A excitation at higher momentum. 

The lower panels correspond to a temperature of T = 100 MeV. The softening of the 
nucleon occupation number due to thermal motion causes a broadening of the three modes 
present in the spectral function. At g = 450 MeV, the ph, Ah and pion peaks are completely 
mixed, although some distinctive strength still prevails at the Ah excitation energy. The 
s-wave and 2p2h self-energy terms completely wash out the structures that were still visible 
at zero temperature. For higher temperatures we have checked that no further structures 
can be resolved in the spectral function, in agreement with Ref. [73]. We however find 
differences at the numerical level, as we expected, since we have implemented different 
hadronic form factors and Landau-Migdal interactions in our model. Moreover, to keep 
closer to phenomenology, we have also implemented in our model the energy-dependent 
p-wave decay width of the A, which favors the mixing of the Ah excitation with the pion 
mode. 

3 Results and Discussion 

3.1 The K meson spectral function in a hot nuclear medium 

We start this section by showing in Fig. [3] the s-wave KN amplitude for 1 = and 1 = 1 
as a function of the center-of-mass energy, Pq, calculated at nuclear matter density, po, and 
temperatures T = MeV (first row) and T = 100 MeV (second row). We have considered 
four different in-medium approaches: a first iteration that only includes Pauli blocking on 
the nucleon intermediate states (dotted lines), and three other self-consistent calculations 
of the K meson self-energy increasing gradually the degree of complexity: one includes the 
dressing of the K meson (long-dashed lines), another considers in addition the mean-field 
binding of the baryons in the various intermediate states (dot-dashed lines) , and, finally, 
the complete model which includes also the pion self-energy (solid lines). Recall that 
the / = amplitude is governed by the behavior of the dynamically generated A(1405) 
resonance. 

We start commenting on the T = results shown in the upper panels, where we clearly 
see that the inclusion of Pauli blocking on the intermediate nucleon states generates the 
/ = A (1405) at higher energies than its position in free space. This has been discussed 
extensively in the literature [3, 5-9] and it is due to the restriction of available phase 
space in the unitarization procedure. The self-consistent incorporation of the attractive 
K self-energy moves the A(1405) back in energy, closer to the free position, while it gets 
diluted due to the opening of the A(1405)A^ —>■ nNA,nNIl decay modes [5,6,8,9]. The 
inclusion of baryon binding has mild effects, lowering slightly the position of the resonance 
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Figure 2: Pion spectral function at p = po for two different momenta and temperatures, 
T = (upper panels) and T = 100 MeV (lower panels). The dashed lines correspond 
to the p-wave self-energy calculation including ph, Ah and short range correlations. The 
solid lines include, in addition, the (small) s-wave self-energy and the 2p2h absorption 
mechanisms. 



peak. A similar smoothing behavior is observed for the / = 1 amplitude as we include 
medium modifications on the intermediate meson-baryon states. As already pointed out in 
Ref. [6,8,9], when pions are dressed new channels are available, such as ANN~^ or T,NN~^ 
(and similarly with AA^^^ components), so the A (1405) gets further diluted. 

At a finite temperature of 100 MeV (lower panels), the A(1405) resonance gets diluted 
and is produced, in general, at lower energies due to the smearing of the Fermi surface that 
reduces the Pauli blocking effects. When pions are dressed the strongly diluted resonance 
moves slightly to higher energies compared to the zero temperature case. The cusp like 
structures that appear at the lower energy side signal the opening of the vrS threshold on 
top of the already opened YNN~^ one. Note that the cusp-like structure appears enhanced 
in the 1 = case. We believe that this is a manifestation at finite temperature and density 
of the two-pole structure of the A(1405). As seen in Refs. [58,60,74-76], this resonance 
is, in fact, the combination of two poles in the complex plane that appear close in energy 
and couple strongly to either KN or ttS states. These two poles move apart and can be 
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Figure 3: Imaginary part of the in-medium KN s-wave amplitude for / = and / = 1 at 
Pq as a function of the center-of-mass energy Pq, for T = 0, 100 MeV, and for the different 
approaches discussed in the text. 



even resolved in a hot medium because density and temperature influence each of them 
differently. Although not shown in the figure, we find that, at twice nuclear matter density, 
the pole that couples more strongly to vrE states moves further below the vrS threshold and 
acquires a clear Breit-Wigner shape. This allows us to conclude that the cusp observed 
in the 1 = amplitude at T = 100 MeV and p = po is essentially a reflection of this 
pole. The shape of the corresponding resonance appears distorted with respect to a usual 
Breit-Wigner because the pole lies just below the threshold of the ttE channel to which it 
couples very strongly, a behavior known as Flatte effect [77]. Note that these structures 
are not seen in the T = results because the ttS threshold in that case is located at around 
1295 MeV, out of the range of the plot. 

Results for the K self-energy and spectral function at po and T = 0, including s- and 
p-wave contributions, obtained in the different self-consistent approaches are compared in 
Fig. m At g = MeV/c, the inclusion of the baryon binding potential (dot-dashed lines) 
moves the quasi-particle peak of the spectral function, which is defined as 



Eqpiqf = q'^ + m\ + Ren^(Egp(g), q\ 



(29) 



to higher energies with respect to the case with no binding (dashed lines). The pion 
dressing (solid lines) further alters the behavior of the self-energy and, hence, that of the 
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Figure 4: Real and imaginary parts of the K self-energy and spectral density, as functions 
of the K energy, for g = MeV/c and q = 450 MeV/c. Results have been obtained at 
T = and po, including (s+p)-wave contributions, for the three self-consistent approaches 
discussed in the text. 



spectral function. The attraction of the antikaon mode decreases, while its width increases 
due to the opening of new decay channels induced by the ph and Ah pion excitations. At 
finite momentum, the p-wave AA^^^, SA^^^ and S*A^^^ excitation modes are clearly visible 
around 300, 400 and 600 MeV, respectively. The latter mode mixes very strongly with the 
quasi-particle peak, making the differences between various self-consistent approaches to 
be less visible in the K spectral function. Similar results were obtained in Ref. [9]. The 
differences between both T = calculations arise mainly by the use in this work of different 
(and more realistic) baryon binding potentials. The effects of the particular details of the 
nucleon spectrum on the K spectral function have been noted recently in the T = study 
of Ref. [78]. 

The effect of finite temperature on the different partial-wave contributions to the K 
self-energy is shown in Fig. |5l In this figure we display the real and imaginary parts of the 
K self-energy together with the K spectral function for the self-consistent approximation 
that dresses baryons and includes the pion self-energy. We show the results as a function 
of the K energy for two different momenta, q = MeV/c (left column) and q = 450 MeV/c 
(right column). The different curves correspond to T = and T = 100 MeV including the 
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Figure 5: Real and imaginary parts of the K self-energy and spectral function at Pq, for 
q = 0,450 MeV/c and T = 0, 100 MeV as function of the K energy, including s- and 
(s + p)-wave contributions in the full self-consistent calculation. 



s-wave and the (s -|- p)-wave contributions. 

According to Eq. ( 125|) . the present model does not give an explicit p-wave contribution 
to the K self-energy at zero momentum. The small differences between the s and the s+p 
calculations observed at g = MeV/c are due to the indirect effects of having included the 
p-wave self-energy in the intermediate meson-baryon loop. The importance of the p-wave 
self-energy is more evident at a finite momentum of g = 450 MeV/c. The effect of the 
subthreshold AA^~^, T,N~^ and E*iV~^ excitations is repulsive at the KN threshold. This 
repulsion together with the strength below threshold induced by the those excitations can 
be easily seen in the spectral function at finite momentum (third row). The quasi-particle 
peak moves to higher energies while the spectral function falls off slowly on the left-hand 
side. 

Temperature results in a softening of the real and imaginary part of the self-energy as 
the Fermi surface is smeared out. The peak of the spectral function moves closer to the 
free position while it extends over a wider range of energies. 

For completeness, we show in Fig. [6] the evolution of the K spectral function with 
increasing temperature for two different densities, po (upper row) and 2po (lower row), 
and two momenta, q = MeV/c (left column) and q = 450 MeV/c (right column), in 
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Figure 6: The K meson spectral function for g = MeV/c and q = 450 MeV/c at po 
and 2po as a function of the K meson energy for different temperatures and for the self- 
consistent calculation including the dressing of baryons and pions. 



the case of the full self-consistent calculation which includes the dressing of baryons and 
the pion self-energy. At g = MeV/c the quasi-particle peak moves to higher energies 
with increasing temperature due to the loss of strength of the attractive effective KN 
interaction. Furthermore, the quasi-particle peak enhances its collisional broadening and 
gets mixed with the strength associated to the A(1405) appearing both to the right (slow 
fall out) and left (cusp-like structures) of the peak. All these effects are less pronounced 
at a finite momentum of g = 450 MeV/c. In this case, the region of the quasi-particle 
peak is exploring K energies of around 700 MeV, where the self-energy has a weaker 
energy dependence with temperature, as seen in Fig. [51 We note that, as opposed to 
the zero momentum case, the width of the quasi-particle peak decreases with increasing 
temperature because of the reduction of the inter-mixing with the S*A^^^ excitations which 
get diluted in a hot medium. As for the density effects, we just note that the quasi-particle 
peaks widens for larger nuclear density due to the enhancement of collision and absorption 
processes. In fact, a significant amount of strength is visible at energy values substantially 
below the quasi-particle peak. The fact that the K spectral function spreads to lower 
energies, even at finite momentum, may have relevant implications on the phenomenology 
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of the meson propagation and decay in a nuclear medium. We will further elaborate on 
this point in the Conclusions section. 

3.2 K meson in nuclear matter at finite temperature 
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Figure 7: The K meson spectral function for g = MeV/c and q = 450 MeV/c at po and 
2po as a function of the K meson energy for different temperatures. 



The evolution with temperature and density of the properties of kaons is also a matter 
of high interest. In particular, kaons are ideal probes to test the high-density phase of 
relativistic heavy-ion collisions at incident energies ranging from 0.6A to 2AGeV and for 
studying the stiffness of the nuclear equation of state [29]. Moreover, there is a strong 
interconnection between the K~^, K~ and channels, which can lead to important changes 
in the 0-meson production in heavy- ion collisions [79]. 

In the S = 1 sector, only the KN channel is available and, hence, the many-body 
dynamics of the K meson in nuclear medium is simplified with respect to the KN case. 
The K spectral function at g = and q = 450 MeV/c is displayed in Fig. [7] for the self- 
consistent calculation that includes the dressing of baryons, for po (upper row) and 2 po 
(lower row) and different temperatures. The K meson is described by a narrow quasi- 
particle peak which dilutes with temperature and density as the phase space for coUisional 
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KN states increases. The s-wave self-energy provides a moderate repulsion at the quasi- 
particle energy, which translates into a shift of the K spectral function to higher energies 
with increasing density. In contrast to the KN case, the inclusion of p-waves has a mild 
effect on the kaon self-energy (compare thin-dashed lines to solid lines at T = 100 MeV and 
q = 450 MeV/c), as they arise from far off-shell YN~^ excitations in crossed kinematics. 
These excitations provide a small, attractive and barely energy dependent contribution to 
the i^ self-energy. 
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Figure 8: The K mass shift for T = MeV as a function of density, obtained within the 
self-consistent calculation and in the Tp approximation. 



We can define the K optical potential in the nuclear medium as 



UK{q.T) 



(30) 



nil 



which, at zero momentum, can be identified as the in-medium shift of the K meson mass. 
The K mass shift obtained in the self-consistent calculation that considers also the nucleon 
binding effects is displayed in Fig. [8] as a function of the nuclear density for T = MeV. 
Our self-consistent results are compared to those of the low-density or Tp approximation, 
obtained by replacing the medium-dependent amplitude by the free-space one in Eq. (^^. 
We observe that the kaon potential at nuclear saturation density in the Tp approximation 
is 4 MeV less repulsive that in the case of the self-consistent approach, which gives a 
repulsion of 29 MeV. This value is in qualitative agreement with other self-consistent 
calculations [32, 80] and close to the 20 MeV of repulsion obtained in K~^ production on 
nuclei by the ANKE experiment of the COSY collaboration [81]. We conclude that the 
low-density theorem for densities below normal nuclear matter is fulfilled within 15 % due 
to the smooth energy dependence of the KN interaction tied to the absence of resonant 
states close to the KN threshold. 
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3.3 In-medium K and K optical potentials at finite temperature 
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Figure 9: The K potential for the full self-consistent calculation at T = 100 MeV and 
0.25po, Po and 2po as a function of momentum. The K potential at T = and po including 
(s -|- p)-waves is also shown. 



In this last subsection we provide the K and K optical potentials at conditions reached 
in heavy-ion collisions for beam energies of the order or less than 2 AGeV, where temper- 
atures can reach values of T = 100 MeV together with densities up to a few times normal 
nuclear density [29,30]. 

In Figures [9] and [TOl we show the K and K optical potentials at T = 100 MeV for 
different densities (0.25po, Po and 2po), including s- (dotted lines) and (s + p)-waves (solid 
lines), as functions of the meson momentum. In the case of po, we also show the potential 
at T = including (s -t-p)-waves (dashed lines). 

The real part of the K potential becomes more attractive as we increase the density, 
going from —4 MeV at 0.25po to —45 MeV at 2po for g = 0, when both s- and p-waves are 
included. The repulsive p-wave contributions to the K potential become larger as density 
increases, reducing substantially the amount of attraction felt by the K. Compared to a 
previous self-consistent calculation using the Jiilich meson-exchange model [7, 8] , here we 
observe a stronger dependence of the optical potential with the K momentum. 
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Figure 10: The K potential for the full self-consistent calculation at T = 100 MeV and 
0.25po; Po S'lid 2po as a function of momentum. The K potential at T = and po including 

(s + p)-waves is also shown. 



The imaginary part of the K optical potential at T = 100 MeV is little affected by p- 
waves, which, as seen in Fig. [5l basically modify the K self-energy below the quasi-particle 
peak. The imaginary part of the potential shows a flat behavior at low momentum and, 
eventually, its magnitude decreases with increasing momentum as the quasi-particle energy 
moves away from the region of YN~^ excitations. 

With respect to the zero temperature case, shown for po in the middle panels of Fig. [HI 
the optical potential at T = 100 MeV shows less structure. The real part amounts to basi- 
cally half the attraction obtained at T = MeV, while the imaginary part gets enhanced at 
low momentum, due to the increase of coUisional width, and reduced at high momentum, 
due to the decoupling of the K quasi-particle mode from the I1*N~^ one. 

The K meson potential changes from 7 MeV at 0.25po to 74 MeV at 2po for g = MeV/c 
and receives its major contribution from the s-wave interaction, the p-wave providing a 
moderately attractive correction. The imaginary part moves from —2 MeV at 0.25po to 
—25 MeV at 2po for q = MeV/c and its magnitude grows moderately with increasing 
momentum as the available phase space also increases. 

Finite temperature affects the real part of the K meson mildly, as can be seen for po 



24 



from comparing the dashed and sohd hnes in the middle panels of Fig. [TOl The magnitude 
of the imaginary part increases with temperature for all momenta, consistently with the 
increase of thermically excited nucleon states. 

From our results for K and K mesons, it is clear that p-wave effects can be neglected 
at subnuclear densities at the level of the quasi-particle properties. However, they become 
substantially important as density increases and are also responsible for a considerable 
amount of the strength at low energies in the spectral function. 

4 Summary, conclusions and outlook 

We have obtained the K and K self-energies in symmetric nuclear matter at finite temper- 
ature from a chiral unitary approach, which incorporates the s- and p-waves of the kaon- 
nucleon interaction. At tree level, the s-wave amplitude is obtained from the Weinberg- 
Tomozawa term of the chiral Lagrangian. Unitarization in coupled channels is imposed 
by solving the Bethe-Salpeter equation with on-shell amplitudes. The model generates 
dynamically the A(1405) resonance in the / = channel. The in-medium solution of the 
s-wave amplitude, which proceeds by a re-evaluation of the meson-baryon loop function, 
accounts for Pauli-blo eking effects, mean-field binding on the nucleons and hyperons via a 
temperature-dependent a — uj model, and the dressing of the pion and kaon through their 
corresponding self-energies. This requires a self-consistent evaluation of the K and K self- 
energies. The p-wave self-energy is accounted for through the corresponding hyperon-hole 
{YN~^) excitations. Finite temperature expressions have been obtained in the Imaginary 
Time Formalism, giving a formal justification of some approximations typically done in the 
literature and, in some cases, improving upon the results of previous works. For instance, 
in this formalism, the Lindhard function of YN~'^ excitations automatically accounts for 
Pauli blocking on the excited hyperons and satisfies the analytical constraints of a retarded 
self-energy. 

The K self-energy and, hence, its spectral function show a strong mixing between the 
quasi-particle peak and the A(1405)A^~^ and YN~'^ excitations. The effect of the p-wave 
YN~'^ subthreshold excitations is repulsive for the K potential, compensating in part the 
attraction provided by the s-wave KN interaction. Temperature softens the p-wave changes 
on the spectral function at the quasi-particle energy. On the other hand, together with 
the s-wave mechanisms, the p-wave self-energy provides a low-energy tail which spreads 
the spectral function considerably, due to the smearing of the Fermi surface for nucleons. 
Similarly, the size of the imaginary part of the potential decreases with momentum, as the 
K mode decouples from subthreshold absorption mechanisms. 

The narrow K spectral function dilutes with density and temperature as the number 
of collisional KN states is increased. A moderate repulsion, coming from the dominant 
s-wave self-energy, moves the quasi-particle peak to higher energies in the hot and dense 
medium. The absence of resonant states close to threshold validates the use of the low- 
density theorem for the K optical potential approximately up to saturation density. The 
inclusion of p-waves has a mild attractive effect on the K self-energy and potential, which 
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results from YN^'^ excitations in crossed kinematics. 

The properties of strange mesons at finite temperature for densities of 2-3 times normal 
nuclear matter have been object of intensive research in the context of relativistic heavy-ion 
collisions at beam energies below 2 AGeV [29] . The comparison between the experimental 
results on production cross sections, energy distributions and polar angle distributions, and 
the different transport-model calculations has lead to several important conclusions, such 
as the coupling between the K~ and the K'^ yields by strangeness exchange and the fact 
that the K~^ and K~ mesons exhibit different freeze-out conditions. However, there is still 
debate on the influence of the kaon-nucleus potential on those observables. The in-medium 
modifications of the K and K properties devised in this paper could be used in transport 
calculations and tested against the data from the current experimental programs in heavy 
ions [29,30]. 

The fact that the K spectral function spreads to low energies, even at finite momentum, 
may have relevant implications on the phenomenology of the (p meson propagation and 
decay in a nuclear medium. The reduced phase space for the dominant decay channel in 
vacuum, -^ KK, makes the meson decay width a sensitive probe of kaon properties 
in a hot and dense medium (the p-wave nature of the (pKK coupling further enhances this 
sensitivity). In [43,44] the meson mass and decay width in nuclear matter were studied 
from a calculation of the K and K self-energies in a chiral unitary framework similar to 
the present work (the most relevant differences and novelties introduced in this work have 
been discussed in previous sections). The overall attraction of the K meson together with 
a sizable broadening of its spectral function (which refiects the fate of the A(1405) in a 
nuclear medium), induced a remarkable increase of F^ of almost one order of magnitude 
at p = po as compared to the width in free space, as several decay mechanisms open in 
the medium such as (pN -^ KY and (pN -^ KttY. The LEPS Collaboration [36] has 
confirmed that the meson width undergoes strong modifications in the medium from the 
study of the inclusive nuclear photoproduction reaction on different nuclei. The observed 
effects even surpass the sizable modifications obtained in [43,44] and predicted for the 
photoproduction reaction in [47]. 

At finite temperature (CERN-SPS, SIS/GSI and FAIR/GSI conditions), despite the K 
peak return towards its free position, we expect a similar or even stronger broadening of 
the meson, as S^ further dilutes in the medium effectively increasing the available phase 
space. Note, in addition, that the presence of thermally excited mesons induces "stimu- 
lated" — > KK decays (as well as diffusion processes) [82,83]. Since Bose enhancement 
is more effective on the lighter modes of the system, the low energy tail of the K spectral 
function may contribute substantially to the decay width. 

At RHIC and LHC conditions, the hot medium is expected to have lower net baryonic 
composition. One may conclude, as a consequence, that the contribution from interactions 
with baryons will be smaller. The relevance of baryonic density effects even at high tem- 
peratures has been stressed for the p and meson clouds in hot and dense matter [2,45,84]. 
In the case of — > KK decays, a finite density of antibaryons allows the K to interact 
with the medium through the charge-conjugated mechanisms described here for the K, 
and viceversa (in the limit of p^ = the K and K self-energies are identical). At small 
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net baryon density, whereas the effective contribution from the real parts of the K, K 
self-energies tend to vanish, it is not the case for the imaginary parts, which are always 
cumulative. Thus even at small baryonic chemical potential, the presence of antibaryons 
makes up for the loss of reactivity from having smaller nuclear densities. Additionally, the 
relevance of kaon interactions with the mesonic gas (ignored in this work) becomes man- 
ifest in this regime, as it has been pointed out in [46,85-87]. This mechanisms, together 
with Bose enhancement of KK decays point towards a sizable increase of the reactivity 
even at very high temperatures. 

Therefore, we plan to study the influence of the KK cloud on the properties of the 
meson in a nuclear medium at finite temperature [88], extending our previous analysis 
for cold nuclear matter [43,44,47]. Such changes on the meson properties are a matter 
of interest in the current and future experimental heavy-ion programs [30,49,50]. In 
particular, the future FAIR facility at GSI will devote special attention to the in-medium 
vector meson spectral functions. The HADES experiment will operate at higher beam 
energy of the order of 8-10 AGeV, providing complementary information on the spectral 
function evolution of vector mesons to the current research program. On the other hand, 
CBM will measure the in-medium spectral functions of short lived vector mesons directly 
by their decay into dilepton pairs. 

With this work we expect to pave the understanding of kaon properties in hot and 
dense matter and provide an essential ingredient for the 0-meson phenomenology in heavy- 
ion collisions. Our results are based on a self-consistent many-body calculation at finite 
temperature which relies on a realistic model of the kaon nucleon interaction thoroughly 
confronted to the kaon nuclear phenomenology. 
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A Finite-temperature Lindhard functions 

We quote here some Lindhard function expressions in the ITF to point out the main 
differences with previous evaluations. 
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A.l YN ^ excitations 

In the ITF, the YN'-^ Lindhard function reads 

J {2iryiLJn + En{p,T)- EY{p + q,T) 

where a;„ is a bosonic Matsubara frequence {iun = VlniiT) and the factor 2 stands for 
spin degeneracy. Consistently with the approximations employed in this work, we have 
only kept the positive energy part of the baryon propagators, while keeping the baryon 
energies fully relativistic and containing also mean-field binding potentials. The nucleon 
and hyperon Fermi distributions, nNy{PiT) = ^e^^'^-^^^''^'^~^'^^^'^ + 1]^"*^, depend on the 
temperature and baryon chemical potential, so that for fixed T and fiB the nucleon and 
hyperon densities are given by 

/d V f d T) 

-—--nN{p,T) , py = uy j -n^'^y^P^^) ' (32) 

with vb the corresponding spin-isospin degeneracy factors, namely z/jv = 4 and z/y = 
(2, 6, 12) for Y = (A, S, S*). All the hyperons are considered as stable particles. 

Note that, differently from the zero temperature case shown below, at finite tempera- 
ture there are occupied hyperon states and therefore the Yh excitation is suppressed by 
hyperon Pauli blocking, as it is evident from the hyperon distribution, uy, which appears 
subtracting in the numerator of the Lindhard function. This can be seen explicitly by 
rewriting Un — Uy as nj^{l + Uy) — '?^y(l + n^). 

Analytical continuation to real energies (iu;„ — > Q'o+i^) gives the following expressions for 
the real and imaginary parts of the finite temperature YN~^ Lindhard function (including 
mean-field binding potentials), 

2tt^ Jo J~i qo + En{p,T) - Ey{p + q,T) 

Imf/yAr-i(go,g;T) = -TT— ^ / dpp 

Zn^ Jo pq 

X [njv(g; T) - nY{p + g, T)]„„ ^(1 - \uo\) 9{qo + En{p, T) - Ey) , 

(33) 

with Uq = UQ{qQ,q,p) = [(go + -E'at — Sy)^ — (My)"^ — p^ — g^]/(2pg), where here g, p refer 
to the modulus of the corresponding three-momentum. 

In the limit of zero temperature and fixed baryon chemical potential, /i^, which then 
coincides with the nucleon Fermi energy, Ep = Ei\i{pp), with pp the Fermi momentum, we 
have 

n^ip, T) -^ riNip) = 9{pp - p) , nY{p, T) ^ . (34) 

The YN~'^ Lindhard function then reads 

UYN-^{qo,q;p) = 2 / 77--^ , ^ x^x ,^ . , 35 

J {2-Kyqo + EN{p)-EY{p + q) + ie 
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and p = 2p^/3tt'^. In [8,43,44] analytical expressions are provided for the non-relativistic 
Fermi gas (i.e., Eb{j)) = Mb + 'p'^ /2Mb above), which we quote here for completeness, 

3 My f ^ ^ 1 



with 



ImUYN-i{(lo,q;p) 

My 

z = 

qpF 



-p < ; 

2 qpF [ 



z')\n 



\z + l\ 



Qo 



4 qpp 



{My - Mn] 



2M 



y 



(36) 



(37) 



We note that terms proportional to {My^ — Mg^) in the denominator of Eq. ( l35l) are 
neglected when doing the angular integration to arrive to Eq. ( l36l) . 

In the literature, one often finds approximate expressions for the finite temperature 
Lindhard function that have been obtained from extensions of the former T = equations. 
In [8] the finite temperature generalization of Eq. (!35|) was obtained by replacing the 
nucleon occupation number, riNij)) = 0{pf —p), by the corresponding Fermi distribution 
at finite temperature, riiqij)) — > nj^{p,T). Analytical expressions (up to a momentum 
integration) read 
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A. 2 ph and Ah excitations 

The evaluation of the ph = NN^'^ Lindhard function at finite temperature and density 
can be found, for instance, in [73,89]. Analytical continuation to real energies from the 
ITF expression leads to 

d^p nN{p, T) - rajv(p + g, T) 



UNN-^{qo,q',T) — un / /^ V? , . .^ . .^ ^ . 

J {2TTyqo + ie + EN{p,T)-EN{p + q,T) 
with un = 4. Similarly, for the Ah Lindhard function we arrive at 



UAN-^{qo,q;T) = ^a J 



d^p 

(2^ 



nN{p,T) -nA{p + q,T) 



qo 



rACgo.q) 



+ Er,{p,T)-EA{p + q,T) 



+ 



nA{p, T) -nN{p + g, T) 



go + i^^i^ + Ea(p, T) - En{p + q, T) 



(41) 
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which we have written exphcitly in terms of direct (AA^~^) phis crossed (A^A^^) contri- 
butions. For convenience, the vrA^A couphng is absorbed in the definition of U/^]^~i and 
thus z/A = ^{fd^/fNY- Note that in Eq. (14T1) we have accounted for the decay width of 
the A resonance. A reahstic treatment of the Ah mechanism should account for the full 
energy dependent A decay width into its dominant channel, nN . We implement the A 
decay width using the formulae in [65], in which Fa only depends on the pion energy and 
momentum (a more detailed study of the in-medium A self-energy at finite temperature 
and density has been reported in [90]). FA(q'o,^) accounts for both direct and crossed 
kinematics and hence the retarded property of f/AAf-i is preserved. We have not accounted 
for binding effects on the nucleon and A in the pion self-energy. First, in the ph excitation 
the baryonic potentials cancel to a large extent. Second, the binding potentials for the A 
resonance are not well known experimentally and we omit them. Therefore, for consistence, 
we do not dress the nucleon in the Ah excitation either. 

The T — >■ limit (at nuclear matter conditions) of Uj^j^-i and U^j^-i can be easily 
obtained with similar prescriptions as in the YN~^ case, namely, n]^{j),T) — > nfq{p) and 
riA -^ 0. Analytic expressions for the non-relativistic ph and Ah Lindhard functions can 
be found in [65]. 

B 5-wave self-energy from Tj^(j^\^ 

We derive in this section a general expression for the contribution to the kaon self-energy 
from the effective in-medium K{K)N scattering amplitude at finite temperature. Let us 
denote by Tj^i^k)n the isospin averaged kaon nucleon scattering amplitude. The kaon self- 
energy, n^(j^)^, follows from closing the nucleon external lines and, following the Feynman 
rules in the ITF, reads 

°° f (Pp 1 

mt^oo-J {27rr iWm - En{p) 

where Un and Wm are bosonic and fermionic Matsubara frequencies, respectively, ia;„ = 
i2mTT and iWm = i(2m + l)7cT + fiB- The sum over the index m is not straightforward 
since Tj^(k)n depends on rri in a non-trivial way. To skip this complication, we can invoke 
a spectral representation for the T-matrix (inherited from the analytical structure of the 
meson-baryon loop function) and we have 

d^p 1 roo ImT^(^)^(f],P;T) 



n^(.).(...,g1T) = -T j: j-^,-jyn 



(27r)3 TV J-oo [iWm " ^7v(p)][ic^n + iWm " ^] 

d'p 1 ^°o,^ ImT^(^)^(^,P;T) ^ 

^^ ~ o , z:. /^N [nN{p, T) - n{n, T)] 



(27r)'^ TT J-oo iuJn — fi + E]\f{p) 

(43) 

with n{ri,T) = [e(^~^s)/'^ + 1]~^ here. The former result, after continuation into the 
real energy axis {iun —* qo + is), provides the thermal kaon self-energy evaluated from 
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the kaon nucleoli scattering amplitude. Note that it includes a Pauli blocking correction 
term, n{Q,T), convoluted with the spectral strength from the imaginary part of the T- 
matrix. At the region in which the principal value of the spectral integration gets its major 
contribution, Q ^ qo + EnIp), the fermion distribution n{Q, T) behaves as a slowly varying 
exponential tail (for the present temperatures under study). We can approximate this term 
by a constant, namely, n{Q,T) ~ n{qQ + En{p),T) and take it out of the integral. The 
dispersion integral over fl then recovers the full amplitude Tf^^^^-^^ and the self-energy can 
be approximated by: 

^k{k)n{% + i^, ?*; r) = y 7^^ Tfi^K)NiQo + En{p),P; T) [n^(p, T) - n{qo + Ejv(p ), T)] . 

(44) 
Note that this procedure is exact for the imaginary part. Eq. fl24l) follows from the former 
result by neglecting the Pauli blocking correction on the fermion degrees of freedom excited 
in the kaon nucleon amplitude (note that in the isospin zero amplitude the strength peaks 
around the A(1405) resonance, and thus one expects this correction to be small with respect 
to that on the nucleon). 
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